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Resumen. Se presenta un metodo numerico para obtener soluciones aproximadas de 
problemas provenientes de la sedimentacion de suspensiones floculadas. Estos procesos 
se usan ampliamente en la industria minera, por ejemplo para recuperar el agua de las 
suspensiones que salen de los procesos de flotacion [ ]. 

La idea principal es aplicar los metodos de multiresolucion a los esquemas desarrollados 
por Burger et al. [2, 3, 4] U observar que el metodo de multiresolucion a describir es 
vital para reducir el costo computacional sin afectar la calidad de la solucion. 



1 INTRODUCCION 

Se introduce el problema fi'sico y su modelacion 
mediante una ley de conservacion fuertemente 
degenerada con flujo no lineal. El efecto de la 
compresibilidad del sedimento puede ser des- 
crito por un termino difusivo fuertemente de- 
generado, mientras el flujo unidimensional con- 
tribuye una discontinuidad de flujo a la ecua- 
cion parcial diferencial. Se presenta un esque- 
ma de segundo orden para resolver este tipo de 
problemas y flnalmente se desarroUan ejemplos 
numericos para comparar con resultados publi- 
cados ([3, 4]). 

En el caso unidimensional, la teon'a de la sedi- 
mentacion produce ecuaciones de equilibrio de 
masa y momentum lineal que pueden simplifl- 
carse hasta obtener una ecuacion de la forma 



sion integrado dado por 



A{u) 



a{s)ds, a(n) ^ 0. (2) 



(1) 



con {x,t) g]0,1[x[0,T[ y el coeficiente de difu- 



Se permite que el coeficiente de difusion a{u) 
sea cero sobre intervalos de u. 
Las soluciones de (1) desarrollan discontinuida- 
des debido a la no linealidad de la funcion de 
densidad de flujo f{u) y a la degeneracion del 
coeflciente de difusion. Esto lleva a consider ar 
soluciones entropicas para tener un problema 
bien puesto. Aun mas, cuando (1) es puramen- 
te hiperbolica, los valores de la solucion se pro- 
pagan sobre rectas caracteristicas que podrian 
intersectar las fronteras del dominio espacio- 
tiempo desde el interior, y esto requiere tratar 
a las condiciones de Dirichlet como condiciones 
entropicas. 

a{u) tiene un comportamiento degenerado, es 
decir, a{u) = para u ^ Uc y a{u) salta en Uc 
a un valor positivo, donde Uc es una constan- 
te llamada concentracion critica. Se enfatiza el 



hecho de que el coeficiente de difusion a{u) es 

degenerado, lo que hace evidente la naturaleza 

hiperbolica-parabolica de (1). 
Considerar el PVIF siguiente 

at" + 9^ (<?(«)" + /(")) = 9^-4(1*), (3:,t) e]0,H[x[0,T[(3) 

u(x,0) = uo(x), xe[0,H], (4) 

u(H,t) = 0, tG]0,T] (5) 

/{u(0,t)) - S:iA(u(0,t)) = 0, tG]0,T], (6) 

conocido como el Prohlema ^. Y el Prohlema 
B 

dtu + a^(q(t)u + f(u)) = 8l^A(u), (x,t) e]o,H{x[o,m) 

u{x,0) = uo{x), xe{0,H], (8) 
q(t)u{H,t) - a^A(u{H,t)) = *(t), tg]0,T] (9) 
/(ti{0, t)) - e^A(u(0, t)) = 0, tg]0,T]. (10) 

Para ambos problemas, / se supone continua 
y diferenciable a trozos, / ^ 0, supp(/) C 
[0,Umax], ||/'||oo ^ oo, a{u) ^ 0, supp(a) C 
supp(/), a(n) = para u ^ Uc, < Uc < Umax, 
q{t) ^ 0, Vi G [0,r], TV{q) < oc, TV{q') < 
oo. En [ ] se prueba la existencia y unicidad 
de solucion entropica para cada uno de estos 
problemas. 

La propiedad de mayor interes, es que general- 
mente se supone el siguiente comportamiento 
para ae{u): 



{ 



— etc., 
> 0, 



U ^ ' 

u > ■ 



{ 



= 0. 
> 0, 



. U > Uc 

(11) 



Este esquema es estable bajo la condicion CFL 
max|/'(u)|^+2max|a(?/)|-^ < 1. (13) 



3 ANALISIS DE MR 

Se presentan los conceptos y definiciones basi- 
cas introducidas por Harten [ ] para el analisis 
de multiresolucion. Considerar el conjunto de 
mallas anidadas diadicas G'', k = 0, . . . , L: 



No 
(14) 



o c> o o o o o o 



o o o 



I 2 3 4 5 6 7 S 9 10 11 12 13 14 15 16 



2 3 



Diferentes escalas de valores puntuales 



4 ALGORITMO DE MR 



1. Inicializacion de pardmetros y variables: 



2 ESQUEMAS DE SEGUNDO ORDEN 

Los terminos advectivo y difusivo son aproxi- 
mados de diferente forma, con el fin de obtener 
una discretizacion que mantenga la conservati- 
vidad. Para la parte advectiva puede utilizarse 
el esquema de Roe clasico con una interpola- 
cion ENO de segundo orden, o bien un esquema 
de Engquist-Osher [7]. Para la parte difusiva, 
se necesita un esquema centrado de segundo 
orden que mantenga la conservatividad. El es- 
quema interior resultante para la ecuacion (1) 
es 



Ax 



Ti+1 n - + 

+ g(nAt)-i+i ^ 

At Ax 

(Aa;)2 



En este caso, se utilizara el 6—limitador 



s; = MM 



con 



Ax ' 2Ax 



Ax 



e e [0, 2], 



mm(a, b, c), si a,b,c > 0, 
MM{a,b,c) := \ max{a,b,c), si a, 6, c < 0, (12) 
0, e.o.c. 



Longitud del dominio H, concentra- 
cion critica Uc, orden de la interpola- 
cion de multiresolucion r, niveles de 
multiresolucion L, mimero de pun- 
tos y paso en la malla fina A'^o y ho, 
y en cada nivel Ni^. y h^, tolerancia 
prescrita e y estrategia de trunca- 
miento e^, constantes de Lipschitz 
para a{u) y f'{u), 

condicion CFL: 

max|/'(n)|- — 1-2 max |a(n)| — y ^ 1. 
u ho u Hq 

paso temporal At, 

CFL ■ ho 



At 



maxu \f'{u) \ +2 max^ | a{u) \ / ho 



estructura de datos. (SPARSE, 
GRADED TREE). 



Aplicacidn de la codificacion a la condi- Se utilizara la funcion cr'g{u) dada por ([4, 5]) 
cidn inicial: Coeficientes de ondelette sig- 



nificativos y posiciones correspondientes. 
Se incluyen los safety points. 



d 

du 



(I00(n/ 



Un 



Evolucion temporal: Se utiliza un metodo Luego 
Runge-Kutta de segundo orden. 



5 EJEMPLOS NUMERICOS 

Se calculan soluciones de los problemas A y B 
utilizando los esquemas numericos descritos en 
la seccion anterior, con una discretizacion para 
el flujo de tipo Enqguist-Osher. Se reproducen 
algunos resultados numericos. 

5.1 SEDIMENTACION BATCH IDEAL 

Proceso de sediment acion batch de suspension 
ideal en una columna de asentamiento: 




l) [Pa\, si n > Uc- 
(16) 

s\ u ^ Uc = 0,23, 
si u > Uc- 

(17) 



* w 


V 


M 


ei 








60 


6.5737 


17.8796 


1.29x10" 


4 


5.33x10" 


h 


1800 (*) 


5.7349 


9.4132 


1.99x10" 


4 


7.42x10" 


5 


3600 (*) 


6.1982 


9.1246 


2.77x10" 


4 


9.61x10" 


5 


7200 (*) 


6.2110 


9.1246 


3.21x10" 


4 


2.41x10" 


4 


14400(*) 


7.9244 


9.4132 


8.92x10" 


4 


6.18x10" 


4 



SuspensioiiGS floculadas, primer ejemplo. Multiresolucion utilizando 



iVo = 129 y i : 









* l=l|i,](MRi 












l = l|h](DF) 
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Izquicrda: Condicion inicial (rayas) y pcrfil dc conccntracion a 
t — 300[s] para cl problcma dc scdimcntacion batch dc suspension 
ideal (Asteriscos). Dcrecha: Coeficientes de ondelette significativos 
correspondientes. 



t[s] 


V 


/J 


ei 




eoo 




60 


4.3457 


7.8456 


2.64x10" 




9.03x10" 


6 


300 


5.6212 


5.8456 


1.70x10" 


5 


1.12x10" 


5 


1800 


5.9443 


14.9168 


7.28x10" 


5 


4.35x10" 




3600 


6.1385 


29.8479 


8.89x10" 


5 


6.50x10" 


5 



Izquierda: Condicion inicial (rayas) y pcrfil de conccntracion a 
t — 3600[s] para el problema de sedimentacion-consolidacion 
(asteriscos). Derecha: Coeficientes de ondelette significativos 
correspondientes. £ — 10""^, Nq — 129 y L — 5. 




Sedimentacion de suspension ideal, e — 1,0 X 10 , A'"q — 257 y L — 5. 

5.2 SEDIMENTACION BATCH CON 
COMPRESION 

Como funcion de densidad de flujo, se utili- 
za una funcion Kynch batch Richardson- Zaki 
con parametros correspondientes a suspension 
de cobre [3]. 

f{u) = -6,05 xW-\{l-u)'^'''[m/s]. (15) 5.3 SEDIMENTACION CONTINUA 



Perfiles de conccntracion hasta t — 12[h], problcma de 
sedimentacion-consolidacion Batch. £ — 10 ^, Nq — 129 y L — 5. 




Pcrfilcs dc conccntracion hasta t — 16[/i] para cl problcma dc 
sedimciitacion continua. s — 5 X 10~^, Nq — 513 y L — 5. 

5.4 REACCION-DIFUSION ID 




Condicion inicial (rayas) y S(u) inicial (puntos), solucion numcrica 
sin multircsolucion (linea), solucion numcrica con multircsolucion 
(asteriscos) y S{u) (puntos-rayas), cn ci ticmpo t — 10 para la cc. dc 
rcaccion-difusion, con ck — 0,8, — 10, L — 7, Nq — 513 ye — 10~^. 
Dctallcs significativos, t — 0,5. 

5.5 TRAFICO VEHICULAR 

MRS solution 








Solucion tridimensional del problcma dc fiujo dc trafico cn una 
rotonda. 

p[cars/mi] MRS solution 

250 n 




Solucion tridimensional del problcma dc flujo dc trafico cn una 
rotonda. 
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